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Abstract 

Beam  misalignments  and  bulk  fluid  velocities  can  influence  the  time  history  and  intensity 
of  laser-induced  thermal  acoustics  (LITA)  signals.  A  closed-form  analytic  expression  for  LITA 
signals  incorporating  these  effects  is  derived,  allowing  the  magnitude  of  beam  misalignment  and 
velocity  to  be  inferred  from  the  signal  shape.  It  is  demonstrated  how  instantaneous,  nonin- 
trusive,  and  remote  measurement  of  sound  speed  (temperature)  and  velocity  (Mach  number) 
can  be  inferred  simultaneously  from  homodyne-detected  LITA  signals.  The  effects  of  different 
forms  of  beam  misalignment  are  explored  experimentally  and  compared  with  theory,  with  good 
agreement,  allowing  the  amount  of  misalignment  to  be  measured  from  the  LITA  signal.  This 
capability  could  be  used  to  correct  experimental  misalignments  and  account  for  the  effects  of 
misalignment  in  other  LITA  measurements.  It  it  shown  that  small  beam  misalignments  have  no 
influence  on  the  accuracy  or  repeatability  of  sound  speed  measurements  with  LITA. 
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1  Introduction 


The  four-wave  mixing  technique  Laser-Induced  Thermal  Acoustics  is  used  in  various  laboratories 
to  measure  the  sound  speed  and  thermal  diffusivity  (e.g.,  Ref.  [1])  as  well  as  fluid  velocities  (e.g.; 
Ref.  [2],  [3]).  A  grating-shaped  pressure  and  temperature  perturbation  is  created  in  a  fluid  by  two 
coherent  intersecting  pulsed  laser  beams.  These  perturbations  evolve  hydrodynamically  in  time  and 
space.  This  evolution  can  be  examined  by  focusing  a  third  continuous  laser  beam  on  the  associated 
density  gratings  at  its  Bragg  angle.  Depending  on  the  modulation  depth  of  the  density  grating,  a 
fraction  of  this  source  beam  is  coherently  scattered  into  a  weak  signal  beam.  The  induced  acoustic 
waves  move  outwards,  modulating  the  density  field  and  hence  the  signal  at  the  grating  Brillouin 
frequency,  or  the  frequency  of  sound  waves  with  the  grating  wavelength.  Hence,  one  approach  to 
measure  the  sound  speed  is  to  measure  the  frequency  of  the  signal  modulation.  However,  finite 
signal  lifetime  and  sampling  resolution  limits  the  accuracy  of  this  technique. 

Another  approach  is  to  model  the  physics  of  the  interactions  that  produce  the  LITA  signal 
and  derive  an  analytic  expression  for  the  signal  as  a  function  of  experimental  parameters,  fluid 
properties  including  the  sound  speed  and  thermal  diffusivity,  and  fluid  velocity.  A  least-squares  fit 
of  this  expression  to  experimental  signals  then  provides  best  estimates  of  the  signal  parameters  and 
hence  fluid  properties.  Cummings  (Ref.  [4])  derived  such  an  expression  that  includes  thermalization 
and  electrostriction  as  the  two  main  mechanisms  for  the  creation  of  the  grating.  However  the 
assumption  was  made  that  the  system  was  in  perfect  optical  alignment  and  that  the  fluid  was  at 
rest.  Obviously  this  approach  breaks  down  when  these  assumptions  are  grossly  violated.  In  these 
cases,  the  theory  is  not  a  valid  representation  of  the  experiment  and  the  fitting  procedure  returns 
erroneous  values. 

The  work  presented  here  extends  the  earlier  theory  to  include  the  effects  of  a  convection  velocity 
and  finite  beam  misalignment.  These  may  then  be  measured  through  the  least-squares  fitting 
procedure,  simplifying  LITA  velocimetry.  Present  velocimetry  techniques  using  LITA  rely  on  the 
measurement  of  a  small  Doppler  shift  of  the  signal  beam  using  heterodyne  detection.  Now,  only 
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from  the  recorded  signal  and  using  the  simpler  homodyne  detection,  sound  speed,  velocities  and 
transport  properties  of  the  fluid  can  be  extracted  simultaneously.  In  addition  we  can  extract  beam 
misalignment  measurements  from  the  signal  and  either  realign  the  optics  with  this  information  or 
take  them  into  account  in  the  data  analysis. 

In  the  first  part  of  this  paper,  the  LITA  analysis  of  Ref.  [4]  is  extended  using  nomenclature 
mostly  adapted  from  that  derivation.  First,  the  grating  creation  due  to  the  opto-acoustic  forcing 
is  modeled.  The  functional  form  of  the  electric  field  grating  in  the  sample  volume  is  derived.  Next, 
the  fluid  response  to  the  electric  field  is  examined.  Since  the  underlying  physics  have  not  changed, 
the  fluid  response  to  electric  fields  has  been  left  unchanged  from  Ref.  [4],  only  the  results  are  given, 
and  the  interested  reader  is  referred  to  there  and  to  Ref.  [5].  Finally,  the  scattering  of  the  source 
beam  into  the  signal  beam  is  modeled  using  the  linearized  equation  of  light  scattering. 

In  the  second  part,  experimental  results  are  presented  that  validate  the  theory. 

2  Frame  of  Reference 

We  define  our  frame  of  reference  so  that  the  origin  of  our  fixed  Cartesian  coordinate  system  is  at 
the  focus  of  the  source  beam  as  shown  in  Fig.  1.  The  source  beam  lies  in  the  x  -  y- plane  and 
forms  an  angle  i p  with  the  x-axis.  The  x  —  2-plane  bisects  the  two  driver  beams  so  that  each  one 
intersects  this  plane  at  an  angle  9.  They  intersect  the  x  -  y- plane  at  angles  cf>i  and  (j>2  (angular 
misalignment).  Their  foci  are  at  r\  and  f*2,  respectively  (spatial  misalignment).  We  denote  the 
time  of  the  short  driver  pulse  as  t  =  0. 

3  Opto- Acoustic  Forcing 

We  first  model  the  electric  field  grating  in  the  sample  volume  defined  by  the  shallow-angle  intersec¬ 
tion  of  the  two  driver  beams.  Assume  each  driver  laser  beam  has  a  Gaussian  profile  with  Gaussian 
half-width  u>  and  denote  the  normalized  electric  fields  of  the  two  driver  lasers  by  E(i\  (?,  t)  and 
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Figure  1:  Frame  of  reference  for  LITA  analysis. 
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Ed  is  the  total  energy  of  the  driver  laser  field,  Pd(t)  is  the  normalized  driver  laser  intensity  history, 
and  Id  is  the  normalized  grating  intensity  distribution.  Terms  neglected  in  Eqn.  2  only  contribute 
a  background  intensity  that  will  have  no  significant  influence  on  the  signal.  The  result  for  Id  is 


Id  = — y  cos  i2kd  [-£  cos  0  +  {y  -  rj)  sin#]} 
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where 


X= 

sin  9 

l  =(xl  +  x2)/2 
Z  =(xi  -x2)/2 


cos  6 

rj={yi  +  y2)l2 

V  =(yi  -  V2)/2 


Z  =  U! 

C  =  (*1  +  z2)/2 
C  =(*1  -Z2)/2. 


(3b) 


At  this  point  we  use  the  fact  the  the  beam  crossing  angle  is  small  and  that  hence  the  X  >  Y,  Z. 
Consequently,  we  neglect  all  variations  in  the  x-direction  including  the  misalignments.  The  side 
effect  of  this  assumption  together  with  the  assumption  of  no  angular  misalignments  is  that  we 
exclude  any  rotation  of  the  grating,  i.e.  it  is  always  perpendicular  to  the  y-direction.  The  proper 
limit  of  Eqn.  3a  is  then 


h  = 


TTLO‘ 


exp{-:f2  [(v  -vf+v2]  -  J2  (z-<0  +(^2 


cos  {q^(y~  r?)} 


(4) 


with  q ^  =  2kd sin#.  Eqn.  4  represents  a  fringe  pattern  perpendicular  to  the  y-direction  with  a 
Gaussian  intensity  profile  centered  at  (0, 77,  (). 

At  a  later  point  we  will  have  to  take  the  3-dimensional  convolution  of  Id  with  the  Green’s 
functions  of  the  fluid  response  to  the  electric  field.  This  operation  will  be  performed  by  multiplying 
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He  for  the  opto-acoustic  response  of  the  fluid  to  thermalization  and  electrostriction  respectively, 
are 


He{q,t)  =Hop\$pi+Hep2$P2  +  Hot$t  +  HqD$d  (6a) 

He{q,t)  =HePi$Pi  +  HeP2§P2  +  HeT$T-.  (6b) 


where 


$F1)2(g,  t)  =  exp  {-Tq2t  ±  icsqt  +  ivtqy  +  iwtqz) 

§T(q,t)  =  exp  {- DTq2t  +  ivtqy  +  iwtqz} 

f)  =  exp  {-(70  +  7ns)t  ~  Dsq2t  +  ivtqy  +  iwtqz}  (6c) 


fi  +  t(A-G)][i-Gn-*n] .  _H * 
Hepi  -  2[i  +  (A-G)2][(i-cn)2  +  n2]  e- 


1 

HeT  =  _[l  +  (A-G)2][l-AIIj 

[1  — 1(7  A  —  G)][l  +  i{  A  —  G)] 
H ^  =  <l  2[1  +  (A  -  G)2] - 


Hod  = 

_  TT* 

—  -tieP2 


n2 

[(i-Gn)2  +  n2][i-An] 

(7  —  1)A 

HeT=  [1  +  (A  -  G)2]  ’ 


(6d) 


and 

n  =  csq/{w  +  1nB  +  Dsq2).  (6e) 

A  is  the  thermal  grating  damping  over  a  period  of  acoustic  motion  and  T  is  the  acoustic  damping 
rate  (r  =  l/2[(7  -  1  )Dr  +  Dv])-  n  is  the  ratio  of  the  sound  wave  frequency  with  the  wave  vector 
qy  to  the  excited  state  energy  decay  rate.  7  is  the  ratio  of  specific  heats,  DT  the  thermal  diffusivity, 
Dv  the  longitudinal  kinematic  diffusivity  of  the  fluid,  and  Ds  the  diffusivity  of  excited-state  target 
molecules.  70  is  the  rate  of  excited-state  energy  decay  caused  by  thermalization  while  7 n9  is  the 
energy  decay  rate  not  due  to  thermalization. 

<f>Pi  and  <I>p9  are  associated  with  the  acoustic  waves  traveling  in  opposite  directions,  is 
related  to  the  density  perturbations  caused  by  the  thermal  grating,  and  <I>D  is  due  to  the  finite 
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thermalization  rate.  A  single-rate  thermalization  process  is  assumed.  The  density-time  behavior  is 
then 

=  - u2IdPd(t )  O  [HgUe  +  Heue)  (7) 

p 

where  Ug  and  He  are  the  approximate  modulation  depths  of  the  thermalization  and  the  electrostric- 
tion  grating  respectively. 

Id(q)  is  peaked  at  q  -  (0,^,0).  So,  approximate  to  first  order,  q  «  qy.  This  implies  the 
assumption  of  plane  wave  fronts  which  is  a  good  approximation  if  the  waves  have  not  traveled  too 
far  from  their  point  of  origin  and  if  the  fringe  spacing  is  small.  Also  note  that  q-  =  0  implies  a 
constant  intensity  in  direct  ion  when  in  fact,  the  intensity  has  a  Gaussian  profile  in  this  direction. 
But  the  error  is  only  a  multiplicative  factor,  since  all  grating  cross  sections  perpendicular  to  qz  are 
self-similar. 

Taking  the  inverse  Fourier  transform  of  the  products 

${P1,P2,T,D}  =  Id{q)§{P\,P2)T,D}{q-t);  (8) 

we  obtain: 


$S,2(r,t)  exp  j-^2  “  ^C2}exp{— Tqlt}cos{(y-  (v±cs)t)q^} 

x  exp  j-^2  [V  ~  (V  +  (u  ±  c«)*)]2  _  £2  [z  ~  +  wt )]  }  (9a) 

=-~  expj-^r?2  -  {~Drqlt}  cos  {{y  -  vt)  qq,} 

x  exp  [y  -  (??  +  (u +  cs)t)]2  - -^2  [2- (( +  ^)]  |  (9b) 

[)(f,t)  =~  exp  |-^772  -  J2C2}  exP  {~Ds<llt  ~  he  +  7 n*)*}  cos  {(y  -  vt)  qit} 

x  exp  [y  -  (y  + (u +  c5)t)]2  - -^2  [2  -  (C  d-^)]  }■  (9c) 

In  Eqns.  9a-c,  the  assumption  is  made  that  the  length  scales  do  not  change  due  to  diffusion  over 
the  LITA  time  scales.  Also,  a  phase  shift  of  rjq ^  has  been  dropped  in  all  three  cosine  terms. 
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5  Acousto-Optical  Scattering 


After  describing  the  creation  and  time  evolution  of  the  density  grating  in  the  sample  volume,  the 
next  step  is  to  model  how  a  continuous  source  beam  focused  at  the  origin  and  incident  on  the 
grating  at  the  phase-matched  angle  ip  scatters  off  the  grating  into  the  coherent  signal  beam. 

If  we  assume  the  density  variations  pi  to  be  small,  we  can  use  the  linearized  equation  of  light 
scattering  (Ref.  [8]).  A  narrow-band  source  beam  E0 (r,  t)  cos  (w0i)  scatters  into  the  electric  field 
Es  by  a  small  disturbance  in  the  susceptibility  x(r,t:u).  In  the  far  field: 

b 2  _  „ 

ES(R,  t;  q)  =  -^^cos(^  '  R  ~  wo^(g,t),  (10) 

where  ks  is  the  wave  vector  of  the  scattered  beam,  u>0  is  the  frequency  of  the  source  beam.  R.  is  the 
position  vector  relative  to  the  scatterer,  and  p{q,t)  is  the  Fourier  transform  of  the  overlap  of  the 
susceptibility  grating  and  the  source  laser  field 

p{r,t)  =  xr(r,t-,u0)Eo{r,t).  (11) 


Just  as  with  the  driver  beams,  assume  the  source  beam  to  have  spatially  a  Gaussian  intensity 
profile.  Hence, 

(12a) 


Es  =  Eo  +  E*0 

Eq  =  ^  exp  {ik0eo  ■  r  -  exp  |  - 


e0  x  r 


a 


(12b) 

(12c) 


eo  =  (cos  ip,  sin  ip,  0). 

Proceeding  as  before  with  the  driver  lasers,  the  electric  field  of  the  source  laser  is  Es{r,t )  = 
Po{t)Io{r)  cos  (ojQt)  with 

nr  (  (  \  2  /  ,  \  2) 

(13a) 


Jo  = 


7r  a* 


exp  s  -  1  r-  I  -  l  — 

U  v 


where  a  is  the  Gaussian  half-width  of  the  source  laser  and 

a 

az  =  <7. 


sin-i p 


(13b) 
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In  Eqns.  13a, b  variations  in  x-direction  have  been  neglected  by  assuming  ox  =  a/ sin-0  ay,az. 
Define  T  D} (D  *0  =  fo(fO${pi  p2rD}>  of  the  overlap  of  the  source  beam  and  the 

evolving  density  grating.  Taking  the  Fourier  transform  of  these  products,  we  can  write 


®{Pl,P2,T,T}(Q'i)  =  •^^r{Pl,P2,r,T}('?:i)S{PnP2,T!T}(^) 


(14a) 


where 


N 


.  1,2  f  1/2  exp  {_  iv  -  X* 


(TUJ 


2  l  r2  +  2a2  }  [Z2  +  2 a2 


y2r 


f  i  yV 

Tpi,2  =exp|~4y2+2ff2  fa-^)2+^ 


2cj? 


T2  +  2a2 


(17  +  (v  ±  cs)t )  (qy  -  <fy) 


x  exp 


1  ZV2 


z  „2 


[  4Z2  +  2a29z+*Z2  +  2a: 

1  y2cr: 


2(7? 


£  (C  +  wt)  q. 


2a,, 


=  exp  ~ 4  y2  +  2j-2  "  9*)  +  *y2  +  2a2  (??  +  (?J/  ~ 


1  Z2a 2  0  .  2cr2 


xexP  -7^2  ,-0^+* 


4Z2  +  2<72^  Z2  +  2<r2 


m  f  p2i  2  (77  +  (v  ±  cs)t)2 

Spi,2  —  exp  <  F^t  y2  +  2cr2  Z2  +  2a2 


C  +  tut )  }>  =  #  £> 

2 


St  —  exp  < 


SD  =  exp  < 


2  2(r/  +  ut)2  2(c  +  tct) 

DtW  y2  +  2a2  Z2  +  2a2 


C  +  wtj  +  iq${v  ±  cs)t 
> 

iqvvt 


2  \  2  (77  +  vt)2  2(C  +  ^) 

—  (70  +  7 nd)i  y2  2CT2  . 


+  2  a2 


4-  iq^vt 


(14b) 


(14c) 


(14d) 


In  Eqn.  14c,  the  contributions  from  a  similar  lobe  centered  at  qy  =  -q^  have  been  neglected.  Using 
this  result  in  Eqn.  10  yields  the  electric  field  of  the  signal  beam. 


Es{q,R,t)  k2ui2 


Po(t) 


AnR  6XP  0  (^s  ’  R  ~  w°  t)}Pd(t) 

o  R  [api^{0)  +  2°]  +  At${t0)  +  AD4>g'0)]  , 


(15a) 
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where  we  used 


Apit2  =UgHgpiy 2  +  UeHep\fi  Ap  —UgHgr  +  UeHep  Ap>  =UgHflD.  (15b) 

6  Detected  LITA  Signal 

The  detected  LITA  signal  using  heterodyne  detection  is  then  simply  the  integral  of  Eqn.  15  over 
the  detection  angle  where  the  detector  is  centered  at  q  =  (0,  <?■*/,,  0).  In  the  limit  for  a  small  detector, 
Eqn.  15  multiplied  by  the  detection  angle  is  the  LITA  signal.  In  the  experiments  presented  here, 
homodyne  detection  is  used.  The  detector  measures  the  intensity  of  the  electric  field,  i.e.  the  square 
of  the  modulus  of  Eqn.  15. 


P2(t)  _16tt2cos2^Ix(  0)1 

X  (AP14di0)  +  AP1^0)  +  AT^0)  +  AD^0)) 

x  +  A*P1^0)*  +  At<S>{t0)*  +  AD${n0)*)  (16) 

In  Eqn.  16  we  used  the  fact  that  the  driver  pulse  («  7  ns)  is  short  compared  to  the  inverse  Brillouin 
frequency  and  it  was  approximated  by  a  Dirac  delta  function  which  eliminated  the  double  temporal 
convolution  over  the  driver  pulse  time  history. 

In  the  large  detector  limit  the  integration  over  the  detection  angle  can  be  approximated  by 
infinite  integrals.  The  final  result  is  then 

I-'hom  kl  |  /  \|2  f  ^2  ^  >2! 

*  exp  0¥)2} [<Pl + mT’ +  D’]  +  (P! + m(T  +  0)1 


Y2  (Y2  +  2a2 


+  exp 


Y2  (Y2  +  2 of) 


Mr  (PiP2*  +  p;p2) 


+  (PiPj*  +  P2P2*  +  TT*  +  TD*  +  T*D  +  DD* 
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where  Pi  =  ApiEpi,  T*  =  etc.,  and  where  the  symbols  Ap\,p2,t.d  and  ^pi,P2,t,d  are 

defined  in  Eqn.  15b  and  Eqn.  14d,  respectively. 

In  the  absence  of  beam  misalignments  and  fluid  velocities.  Eqn.  17  collapses  to  the  solution 
in  Ref.  [4]  with  the  exception  of  a  multiplicative  constant.  This  difference  is  due  to  the  fact  that 
variations  in  the  ^-direction  were  dropped  much  earlier  than  in  Ref.  [4]. 

In  the  absence  of  bulk  fluid  velocities,  only  the  misalignment  rj  changes  the  time  history  of  the 
LITA  signal  whereas  77,  £,  and  £  simply  decrease  the  signal  intensity.  Also,  Eqn.  17  is  symmetric 
w.r.t.  all  possible  misalignments.  In  the  absence  of  beam  misalignments,  both  fluid  velocities  have 
a  similar  influence  on  the  signal  that  resembles  additional  acoustic  damping  and  thermal  diffusion. 
Among  all  misalignment  and  velocity  components,  77  and  v  are  the  most  interesting.  In  particular 
by  introducing  a  known  rj  it  is  possible  to  measure  the  time  it  takes  the  density  grating  to  travel  to 
the  interrogation  beam,  thereby  obtaining  the  velocity  component  v  more  precisely  than  by  means 
of  an  enhanced  decay  rate.  Fig.  2  shows  how  intentionally  misaligning  the  beans  can  produce  a 
dependence  of  the  signal  on  velocity  that  is  more  orthogonal  to  other  parameters  and  therefore  more 
accurately  inferred  from  signals.  Note  that  oscillations  are  still  present  in  the  signal  so  that  the 
sound  speed  can  also  be  measured.  Depending  on  the  range  of  velocities  to  be  measured,  suitable 
values  for  77,  a;,  and  9  can  be  chosen. 

7  Experimental  Setup 

Figure  3  depicts  a  schematic  diagram  of  the  experimental  setup.  A  Q-switched,  frequency-doubled 
Nd:YAG  laser  (Spectra  Physics,  GCR-150-10)  drives  a  dye  laser  emitting  ~  10-mJ,  7-ns  pulses 
at  589  nm  with  a  repetition  rate  of  10  Hz.  Behind  iris  i2,  beam  splitter  bs  splits  the  beam  in 
approximately  equal-intensity  halves.  Beam  splitter  bs  and  mirrors  m3,  m4,  m5  are  placed  to 
match  the  path  length  of  the  two  beams  within  ~  1  mm.  Lens  12  with  focal  length  of  750  mm 
focuses  both  beams  onto  the  sample  volume  in  the  test  section.  Lens  12  and  mirror  m6  are  mounted 
on  a  translation  stage  with  a  10-^m-resolution  micrometer  drive.  The  test  section  is  a  high-pressure 
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Figure  2:  Theoretical  LITA  signals  for  atmospheric  air  from  Eqn.  17  with  0  =  1.23°,  uj  =  370  fim, 
a  =  700  fim,  misalignment  rj  =  —  2ui  and  with  in  fluid  flow  in  y-direction  with  different  Mach 
numbers;  a)  M  =  0,  b)  M  =  0.25,  c)  M  =  0.5,  d)  M  =  0.75,  e)  M  =  1.0,  f)  M  =  1.5,  g)  M  —  2.0, 
h)  M  =  3.0. 
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Figure  3:  Schematic  diagram  of  experimental  LITA  setup. 

bomb  with  optical  access  via  anti-reflective  coated  BK7  windows  on  opposite  sides. 

A  CW  Argon-ion  laser  (Spectra  Physics,  Model  165)  at  488  nm  provides  0.5  W  for  the  source 
beam.  A  chopper  wheel  ch  (Scitech  Instruments,  Optical  Chopper)  blocks  the  beam  except  for 
pulses  of  1-ms  duration  which  are  synchronized  (LabSmith  LC880)  with  the  Q-switch  trigger  of 
the  Nd:YAG  laser.  As  with  the  driver  laser,  iris  il  partially  removes  unwanted  beam  modes,  and 
lens  11  (f=l  m)  focuses  the  beam  onto  the  sample  volume.  The  beam  passes  just  over  mirror  m6 
which  directs  the  driver  beams  into  the  test  section.  Note  that  this  optical  setup  introduces  a  small 
angular  beam  misalignment  with  4>\  —  fo-  Mirror  m7  directs  the  scattered  beam  through  iris  i3 
into  the  receiver  unit  where  it  is  focused  on  and  passed  through  pinhole  ph  with  a  diameter  of  400 
pm. 

The  signal  beam  is  detected  by  a  photomultiplier  tube  (Hainmamatsu  model  OPTO-8)  and 
recorded  on  a  digital  storage  oscilloscope  (Tektronix,  TDS  640A)  from  which  it  is  transferred  over 
an  IEEE  488  bus  to  a  personal  computer  for  data  analysis.  The  transmitted  part  of  the  source 
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(interrogation)  beam  is  blocked  as  is  one  of  the  two  driver  beams  exiting  the  test  section.  The  other 
is  detected  by  photo  detector  PD  (Thor  Labs,  DET-2SI)  and  used  to  trigger  the  data  acquisition. 
Each  signal  contains  2,000  data  points  taken  at  a  sampling  rate  of  500  MS/s. 

8  Procedure 

For  all  experiments  presented  here,  the  test  section  was  filled  with  atmospheric  air  seeded  with 
N02  at  concentrations  on  the  order  of  parts  per  million.  The  low  level  of  seeding  enhances  the 
signal  level  by  typically  two  orders  of  magnitude  without  changing  the  results  of  the  measurements. 
Also,  this  seeding  makes  thermalization  predominant  over  electrostriction  which  can  therefore  be 
ignored. 

For  the  data  analysis,  a  standard  personal  computer  (Pentium,  150MHz)  is  used  to  find  the 
least-squares  fit  of  a  theoretical  signal  (Eqn.  17)  to  an  experimental  trace.  To  keep  the  number  of 
floating  fitting  parameters  as  low  as  possible,  values  for  the  beam  diameters  to  and  a,  beam  crossing 
angle  9,  etc.  are  obtained  from  a  set  of  calibration  measurements.  A  Levenberg-Marquart  scheme 
is  used  for  the  least-squares  fit,  which  is  a  combination  of  the  inverse  Hessian  (multidimensional 
form  of  Newton’s  method)  scheme  and  the  method  of  steepest  descent  (Ref.  [9]). 

Eqn.  17'  was  validated  by  moving  the  translation  stage  with  lens  12  and  mirror  m6  in  the  y- 
direction,  thus  creating  a  known  misalignment  rj.  The  measurement  started  at  a  value  of  rj  =  -2 
mm  and  a  trace  was  recorded  every  10  //m  until  reaching  rj  =  +2  mm.  Every  trace  was  averaged 
over  64  driver-laser  shots  to  reduce  the  noise  levels  at  larger  misalignments. 

For  these  measurements  in  which  the  correct  beam  misalignment  was  to  be  inferred  from  the 
signal  shape,  only  the  misalignment  component  in  question  and  Ug,  the  thermal  grating  modulation 
depth,  were  adjusted  during  the  numerical  fit.  The  latter  parameter  had  to  be  included  to  serve  as 
a  multiplicative  factor  since  Eqn.  17  cannot  accurately  give  the  absolute  signal  but  rather  a  relative 
time  history  whose  total  amplitude  depends  additionally  on  the  characteristics  of  the  detector  and 
other  factors. 
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In  a  second  set  of  measurements,  we  wanted  to  find  whether  small  beam  misalignments  affect  the 
repeatability  and/or  accuracy  of  LITA  measurements.  This  was  of  particular  interest  since  it  did 
not  seem  not  possible  to  detect  such  misalignments  by  the  fitting  procedure.  Hence,  at  a  number  of 
fixed  beam  misalignments  77,  500  64-shot  averages  were  obtained  as  data  for  a  statistical  analysis. 
Here,  only  the  sound  speed,  thermal  diffusivity,  and  Uq  where  adjusted  during  the  fitting  procedure. 
These  data  sets  were  analyzed  three  different  ways:  first,  with  all  misalignments  (wrongly)  set 
to  zero,  secondly  with  the  misalignments  held  fixed  at  their  correct  values,  and  finally  with  the 
misalignment  rj  set  as  a  fitting  parameter  during  the  data  analysis. 

9  Results,  Discussion 

Fig.  4  shows  LITA  signals  recorded  with  various  misalignments.  For  large  values  of  rj  (Fig.  4a  &  b), 
the  only  visible  signal  is  due  to  an  isolated  acoustic  wave  passing  through  the  source  beam.  We  see 
that  the  onset  of  the  signal  has  a  steeper  slope  than  the  tail,  a  result  of  acoustic  damping.  Beam 
diameters  can  be  inferred  accurately  from  these  traces.  The  width  of  the  hump  in  the  signal  gives 
a  good  measure  of  the  driver  beam  width  u.  In  the  absence  of  acoustic  damping,  the  deviation 
of  the  exact  shape  of  the  signal  from  a  Gaussian  profile  gives  information  about  the  source  beam 
diameter  a.  But,  since  F  is  known  in  these  experiments,  this  information  can  still  be  determined 
using  the  least-squares  fit  to  Eqn.  17. 

In  Fig.  4c  we  see  the  first  oscillations  in  the  signal.  This  indicates  that  at  this  location  there  is 
a  sizable  overlap  of  the  source  beam  with  both  the  thermal  grating  and  the  acoustic  waves.  As  rj 
is  decreased  further,  the  hump  moves  to  earlier  times  and  the  oscillations  become  stronger.  Note 
the  subtle  change  in  signal  shape  between  Figs.  4e  &  f  considering  the  large  step  in  rj. 

Fig.  5  shows  an  experimental  LITA  signal  together  with  the  result  of  the  fitting  procedure. 
For  clarity,  they  are  plotted  separately.  We  see  how  precise  Eqn.  17  captures  the  features  of  a 
misaligned  signal. 

Fig.  6  shows  the  beam  misalignment  inferred  by  the  fitting  routine  vs,  the  actual  misalignment 
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Figure  4:  LITA  signals  of  atmospheric  air,  a  =  700  ^m,  u  =  370  jum.  6  =  1.23°,  and:  a)  77 
1100  /xm,  b)  rj  =  900  ^m,  c)  fj  =  700  gm,  d)  fj  —  600  jtim,  e)  fj  =  500  /im,  f)  fj  —  0  fjm. 


Figure  5:  Experimental  signal  (top)  and  fitted  theoretical  signal  (bottom)  for  ry  =  500  //m  with 
iv  =  370  urn,  a  ~  700  um,  6  =  1.23°  in  atmospheric  air. 


set  by  the  translation  stage.  The  dashed  lines  with  a  slope  of  ±1  correspond  to  the  ideal  case.  For 
rj  >  550  /i m,  the  extracted  values  for  rj  had  the  correct  slope  but  were  generally  too  large  by  about 
200  /im.  For  very  large  values  of  rj,  the  larger  scatter  in  the  measurements  is  due  to  the  increased 
noise  level  in  the  signals.  A  linear  regression  for  the  region  \rj\  >  600  /mi  gives  slopes  in  these 
regions  of  to  =  -  .946  and  to  =  0.995,  respectively,  with  a  coefficient  of  determination,r2,  greater 
than  0.98.  The  standard  error  for  the  measured  rj  is  22  /dra¬ 
in  the  center  region,  the  detected  misalignment  forms  a  plateau  at  300  /zm.  The  reason  for 
this  behavior  is  probably  that  the  intensity  profiles  of  the  beams  in  the  experiments  are  not  well 
represented  by  Gaussians.  Fig.  7,  which  shows  the  peak  signal  intensity  vs.  fj,  also  provides  evidence 
that  this  assumption  is  not  satisfied.  If  we  give  the  fitting  routine  the  option  of  an  additional  fitting 
parameter  whose  effect  on  the  signal  depends  strongly  on  the  beam  profiles,  e.g.,  misalignment, 
while  not  accurately  modeling  the  beam  profiles,  it  is  unlikely  that  the  fitting  routine  will  return 
the  proper  value  of  the  parameter. 

The  dashed  line  in  Fig.  7  represents  the  theoretical  value  of  the  (visible)  peak  intensity.  The 
intensities  are  normalized  with  the  intensity  at  f?  =  0.  First,  consider  the  theoretical  curve  in  which 
we  can  distinguish  three  regimes. 

For  \fj\  <  650  /im  (regime  I)  we  find  a  Gaussian  behavior  for  the  signal  intensity  w.r.t.  rj.  Here, 
the  peak  value  occurs  early  on  in  the  LITA  signal  (Figs.  4e  &  f)  due  to  constructive  interference 
of  the  thermal  grating  with  the  acoustic  waves.  Both  gratings  have  a  Gaussian  profile  so  the  peak 
intensity  scales  the  same  way. 

In  regime  II  (650  /im  <  |7?|  <  1300  /im)  the  peak  intensity  occurs  at  later  times  in  the  signal 
when  an  isolated  acoustic  wave  passes  the  interrogation  beam  (Figs.  4a  &  b).  The  travel  time 
of  the  wave  scales  linearly  with  the  misalignment.  The  wave’s  amplitude  decreases  exponentially 
with  the  acoustic  damping  rate  V.  This  explains  the  linear  behavior  on  a  logarithmic  plot  for  large 
misalignments  in  Fig.  7.  Figs.  4c  &:  d  mark  the  transition  between  these  two  regimes. 

For  \fj\  >  1300  /im  (regime  III)  the  theoretical  signal  intensity  decreases  again  faster  than  in 
region  II.  This  is  an  artifact  of  the  limited  duration  of  the  data  acquisition.  In  this  regime,  the 
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Figure  6:  Measured  vs.  true  misalignment  rj. 
The  dashed  lines  have  a  slope  of  ±1  respec¬ 
tively.  The  experimental  conditions  are  the 
same  as  in  Fig.  4. 


true  <r|>  [jim] 

Figure  7:  Normalized  peak  signal  intensity  vs. 
rj.  The  experimental  conditions  are  the  same 
as  in  Fig.  4. 


peak  of  the  hump  visible  in  Figs.  4a  &  b  occurs  after  the  end  of  the  recording  time. 

The  measured  peak  intensities  are  higher  than  predicted  over  a  range  of  misalignments  near 
perfect  alignment  before  they  approach  the  theoretical  values  at  larger  rj.  The  most  likely  explana¬ 
tion  is  the  presence  of  laser  modes  other  than  TEM0o-  We  can  see  the  transition  from  regime  I  to 
regime  II  in  Fig.  7.  The  transition  to  III  is  less  pronounced  as  at  this  point  signals  are  very  weak. 
We  can  see  the  effect  of  noise  from  the  large  amount  of  scatter  in  the  intensities  in  Fig.  7. 

As  evidenced  by  Fig.  6  it  was  not  possible  to  extract  accurate  measurements  of  small  misalign¬ 
ments  from  the  signal  shape.  This  finding  prompts  the  question  of  whether  small  misalignments 
have  an  impact  on  the  accuracy  and  uncertainty  of  other  LITA  measurements.  Figs.  8a  &  b  re¬ 
spectively  show  the  uncertainty  and  the  accuracy  of  LITA  measurements  of  sound  speed  vs.  rj  for 
the  three  different  fitting  strategies. 

All  three  fitting  strategies  yield  almost  identical  results.  For  \rj\  <  300  /jm  the  uncertainty 
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Figure  8:  Uncertainty  and  error  of  sound  speed  vs.  tj.  The  sound  speed  for  77  —  0  is  taken  as 
reference.  Three  different  fitting  strategies  were  used.  In  the  strategy  denoted  by  “float,”  rj  is  a 
floating  fitting  parameter,  in  “fixed  0.'’  fj  is  fixed  at  zero,  and  in  fixed  corr,  t]  was  held  constant 
at  the  correct  preset  value.  The  experimental  conditions  are  the  same  as  in  Fig.  4. 
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remains  approximately  constant  at  0.!%  as  does  the  error  of  0.2%,  where  the  measured  sound 
speed  from  the  data  set  with  Tj  =  0  has  been  taken  as  reference  value.  Only  for  M  >  300  pm  do 

the  uncertainty  and  error  increase  significantly. 

Other  types  of  misalignments  that  have  been  carried  ou,  experimentally  include  measurements 

where  only  one  driver  beam  was  displaced  in  y-direction,  thus  inducing  a  simnltaneous  misalignment 
in  v  and  moving  both  beams  in  opposite  directions  (only  V)  and  measurements  with  misalign¬ 
ments  in  the  z-direction,  i.e.,  sweeps  through  <  and  <.  Misalignments  in  the  ^-direction  produced 
similar  results  to  those  explained  in  detail  above.  For  nonzero  values  of  <  or  ?,  no  change  in  signal 
shape  has  been  observed  as  Eqn.  17  predicts. 

When  a  small-diameter,  40  pm  pinhole  ph  was  used  in  the  receiver  unit,  another  form  of 
misalignment  was  observed  that  is  not  included  in  Eqn.  17.  Changes  in  signal  shape  did  not  appear 
to  be  symmetric  w.r.t.  y.  The  most  likely  explanation  for  this  asymetry  is  that  part  of  the  s.gnal 
beam  was  occluded  by  the  pinhole.  An  illustration  of  this  effect  is  shown  in  Fig.  9.  In  the  upper 
half  we  see  the  intensity  of  the  overlap  of  the  density  grating  with  the  source  beam  in  physical  space 
along  the  y-axis  as  it  develops  over  time.  The  bright  and  dark  stripes  correspond  to  the  oscillations 
in  the  early  stages  of  the  signal  whereas  the  bright  tail  reflects  the  exponentially  decaying  tad  of 

the  LITA  signal. 

In  the  large  detector  limit  we  assume  that  the  detector  integrates  over  the  entire  horizontal 
range  in  Fig.  9  (and  beyond).  Imagine  that  half  of  the  signal  beam  is  blocked  so  the  receiver  sees 
only  the  left  or  right  half  of  the  signal  beam  (receiver  misalignment).  The  resnlting  LITA  signals 
are  given  in  the  lower  half  of  Fig.  9.  We  see  that  for  1,  =  0  this  has  no  influence  on  the  signal  shape. 

In  fact,  we  could  take  any  region  along  the  y-axis  and  would  end  up  with  the  same  signal  shape, 
meaning  that  this  kind  of  misalignment  cannot  be  detected  if  the  optics  are  otherwise  perfectly 

aligned. 

However,  if  V  *  O.then  the  time  histories  at  different  positions  y  are  no  longer  self  similar, 
resulting  in  different  signals  depending  on  which  part  of  the  signal  beam  the  detector  sees. 

Figure  10  provides  evidence  that  detector  misalignment  produces  asymmetry  in  beam-misalignment 
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Figure  9:  Intensity  of  overlap  of  density  grating  with  source  beam  from  the  inverse  Fourier  Trans¬ 
form  of  Eqn.  14  for  atmospheric  air  and  at  z  =  0,  i.ev  along  the  grating  center  line.  On  the  top, 
the  time  history  as  a  function  of  y-position  is  plotted  -  on  the  left  for  perfect  beam  alignment,  on 
the  right  with  fj  =  u.  The  experimental  parameters  are:  a  =  u>  =  500  pm,  A 4  =  589  nm,  9  =  0.8°. 
White  corresponds  to  high  intensities,  black  to  low  intensities.  In  the  lower  portion,  the  resulting 
LITA  signals  are  shown  if  the  receiver  only  detects  the  range  y  <  0  or  only  y  >  0,  respectively. 
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Figure  10:  Two  traces  from  a  sweep  through  values  of  fj.  Left  trace:  fj  =  -650  p m;  right  trace: 
fj  —  +650  p m.  A  small  pinhole  ph  with  a  diameter  of  only  40  pm  in  the  receiver  unit  blocked 
portions  of  the  signal  beam.  Compare  with  theoretical  results  in  Fig.  9. 

effects.  Here,  two  traces  taken  at  rj  —  ±650  pm  are  shown  from  earlier  measurements  with  the 
40  pm  pinhole.  We  clearly  see  the  difference  in  the  signal  shape.  Compare  Fig.  10  with  the  two 
traces  in  the  lower  right  of  Fig.  9. 

To  fully  verify  the  origin  of  the  asymmetry,  we  performed  a  series  of  measurements  using  the  40- 
pm  pinhole  in  which,  for  several  values  of  fj.  the  mirror  ml  was  translated  to  block  different  regions 
of  the  signal  beam.  The  influence  of  receiver  misalignment  on  the  signal  shape  was  demonstrated, 
but  since  this  form  of  misalignment  is  not  included  in  Fig.  17,  no  quantitative  data  can  be  given. 
Subsequently,  the  use  of  a  larger  size  pinhole  eliminated  the  asymmetry  at  the  cost  of  increased 
sensitivity  to  incoherently  scattered  light  and  luminosity. 

Finally,  consider  the  case  of  angular  beam  misalignments  and  misalignments  in  the  x-direction. 
Both  can  produce  a  rotation  of  the  principal  axes  of  the  grating  structure.  As  a  consequence,  the 
phase  matching  condition  for  the  source  beam  may  no  longer  be  satisfied  resulting  in  a  weaker 
signal  beam  that  may  even  be  scattered  into  a  different  direction  than  anticipated.  This  alone, 


25 


-  V  0 


however,  does  not  change  the  shape  of  the  signal  as  long  as  the  entire  signal  beam  falls  on  the 
detector. 

10  Conclusions 

An  analytical  expression  for  the  magnitude  and  time  history  of  LITA  signals  from  finite  Gaussian 
beams  in  the  presence  of  fluid  velocities  and  the  most  common  forms  of  beam  misalignment  has 
been  derived.  In  experiments,  some  deviations  from  this  expression  were  observed  that  were  due 
to  non-Gaussian  laser  beam  profiles.  In  regions  where  this  assumption  was  satisfied  however,  the 
experiments  showed  very  good  agreement  with  the  theory. 

It  was  demonstrated  how  beam  misalignment  can  be  detected  quantitatively  only  from  the 
shape  of  the  signal.  Small  misalignments  that  could  not  be  accurately  measured  were  shown  not  to 
influence  the  accuracy  or  repeatability  of  sound  speed  measurements,  even  when  the  data  analysis 
incorrectly  assumed  that  no  misalignments  were  present  during  the  measurement. 

Controlled  beam  misalignments  may  be  advantageous  for  optimizing  the  accuracy  of  measure¬ 
ments  of  a  particular  parameter.  For  example,  beam  misalignments  can  improve  the  accuracy  of 
homodyne  LITA  velocimetry. 
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